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Abstract. We present the convergence analysis of the cell average technique, introduced in 
|12j , to solve the nonlinear continuous Smoluchowski coagulation equation. It is shown that the 
technique is second order accurate on uniform grids and first order accurate on non-uniform 
smooth (geometric) grids. As an essential ingredient, the consistency of the technique is thor- 
j^jQj oughly discussed. 
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1. Introduction 

In this article we study some mathematical issues related to the convergence of the cell average 

technique (CAT) for solving the continuous Smoluchowski coagulation equation (SCE) which 

describes the dynamic evolution of particle growth. This model has many applications in biology, 

polymer science, astrophysics and oil industry etc. The nonlinear continuous SCE reads as 

df(t r) 1 f x f°° 
(1.1) -1±L-L = -J P(x-y,y)f(t,x-y)f(t,y)dy- J (3(x,y)f(t,x)f(t,y)dy, 

with 

f(x,0) = f in (x) > 0, x€]0,oo[' 

Here the number density of particles of volume x > at time t > is denoted by f(x,t). 
The coagulation kernel /3(x,y) represents the rate at which particles of volume x coalesce with 
particles of volume y. It will be assumed throughout the article that /3(x,y) = f3(y,x) for all 
x, y > 0, i.e. symmetric and j3{x, y) = for either x = or y = 0. The integrals on the right-hand 
side of (jl.ip represent, respectively, 

• birth of particles of volume x as a result of coagulation events of particles with volumes 
y and x — y (0 < y < x) 

• death of particles of volume x due to the coagulation events with particles of volume y 
(0 < y < oo). 

There are several results available on the existence and uniqueness of solutions to (jl.ip . see e.g. 
[D O HI O El HOI HH [IS HSl HH UHl [20]. To show all these results, one always needs certain 
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growth conditions on the kernels. The SCE (jl.ip is analytically solvable only for some specific 
examples of kernels, see [21 E]. Because of these restrictions, we are always interested to 
develop new numerical techniques with a detailed study of their mathematical analysis. Among 
all numerical methods for solving SCE (jl.ip . the sectional methods are widely used because 
they are computationally attractive. Recently, the cell average technique is introduced in |12] 
which preserves all advantages of the existing sectional methods. Unlike the previous sectional 
methods, it gives very accurate prediction of selected higher moments and also provides quite 
satisfactory numerical results for number density. This gives us a strong motivation to analyze 
this technique mathematically. To apply a numerical method, first we need to consider the 
following truncated form of the problem (jl.ip by taking a finite computational domain ]0,R] 
where < R < oo. 

dn(t,x) 1 f x f R 

(i.2) — — = 2j ^ X ~ ~ ~ J P(x,y)n(t,x)n(t,y)dy, 

with 

n(x,0) = n m (x) > 0, x€Sl:=]0,R], 

where n(t, x) represents the solution to the truncated equation (|1.2p . The existence and unique- 
ness of non-negative solutions for the truncated SCE (|1.2p has been shown in [21 120] . In 
|U IH HI [101 [151 [SO], it is proven that the sequence of solutions to the truncated problems 
converge weakly to the solution of the original problem in a weighted L l space as R — > oo for 
certain classes of kernels. 

The purpose of this work is to demonstrate the convergence analysis of the cell average 
technique for solving SCE (|1.2p on uniform and non-uniform smooth geometric grids. To the 
best of author's knowledge, this is the first attempt to show the convergence of CAT for solving 
nonlinear continuous SCE. The work presented here is motivated from [13] where the convergence 
of CAT is introduced for solving linear continuous breakage (fragmentation) problems. 

The plan of this paper is as follows. The mathematical formulation of CAT is recalled in the 
next section. At the end of Section [21 the main convergence theorem 12.41 is proved. To fulfill 
the requirements of main result, the consistency of the method and Lipschitz conditions are 
investigated in Section [3] and HI respectively. Finally, some conclusions are made in Section [SJ 

2. The cell average technique 

The cell average technique approximates the total number of particles in finite number of 
cells. As a first step, the continuous interval £1 :=]0, R] is divided into a small number of cells 
defining size classes 

Aj ■■=}x i ^i/2,x i+1 / 2 }, i = l,...,I, 

with 

Xl/2 = 0, X/+1/2 = R- 
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The representative of each size class, usually the center point of each cell Xj = (xj_ 1 / 2 +^i+i/2)/2, 
is called pivot or grid point. We introduce Ax m i n and Ax to satisfy 

Ax min < Axi = x i+1/2 - aJi-i/2 < Ax. 

For the purpose of later analysis we assume quasi uniformity of the grids, i.e. 

*- lJ 'inm 

where K is a positive constant. The total number of particles in the ilh. cell is given as 

fXi+l/2 

(2.2) Ni(t) = / n(t,x)dx. 

Integrating the continuous equation (|1.2p over the ith cell we obtain 

^ = B,-D„ i = l,...,I. 
at 

The total birth rate B{ and the death rate Di are given as 

(2.3) Bi = - / (3(x-y,y)n(t,x-y)n(t,y)dydx, 

z Jxi-i/2 Jo 

and 

r x i+\/2 

(2.4) Di= / p(x,y)n(t,y)n(t,x)dydx. 

The above equations yield a semi-discrete system in M. 1 

dN 

(2.5) — = B - D, with N(0) = N m , 

at 

where N,B,D € M 7 . The ith component of vectors N,B, and D are respectively defined in 
(|2.2p - (|2.4p . The vector N is formed by the vector of values of the step function obtained by 
I? projection of the exact solution n into the space of step functions, which are constant on 
each cell. Note that this projection error can easily be shown of second order, see [7j. The total 
discrete birth and death rates of particles are evaluated by substituting the number density 
approximation 

n(t,x) w5^JVi(t)5(x-Xi) 



i=i 



into equations (|2.3|) and (|2.4|) as 

£ 

Xi-l/2< x j+ x k<Xz + l/2 



(2.6) Bi= U-lkkjP( x k^M(t)N k (t), 



and 

I 

(2.7) D i = N i (t)^2(3(x i ,x j )N J (t). 
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Here Bi and L>i denote the discrete birth and death rates, respectively, in the ith cell. The total 
volume flux Vi into cell i as a result of aggregation is given by 



I f x i+l/2 f x 

(2.8) v i = o / x ^ x - y 

Z JXi-1/2 J 



,y)n(t,x - y)n(t,y)dydx. 



Similarly to the discrete birth rate the discrete volume flux can be obtained as 

j>k 

(2.9) Vi= (l-^Sj^Pix^x^N^N^ixj + Xk). 

x i-l/2< x j+ x k< x i+l/2 

Consequently, the average volume Vi € [aJi-i/2> #1+1/2] °f au new born particles in the ith cell 
can be evaluated as 

(2.10) Vi = ^, Bi>0. 

Bi 

We do not need volume average vi in case of Bi = 0. However, for Bi = 0, we can fix Vi = Xj. 
Here we consider that all new born particles in the ith cell are assigned temporarily to the 
average volume Vi. If the average volume vi is same as the pivot point xi then the total birth Bi 
of the new born particles can be assigned to the pivot x\ only. But this is rarely possible, and 
hence, the total birth Bi has to be assigned to the neighboring pivots in such a way that the 
total number and mass remain conserved during this reassignment. Finally, the resultant set of 
ODEs takes the following form 

(2.11) d -§ = B? A -b? A . 

The above discretized system can also be written in the following vector form 

dN ...... 

(2.12) — = B(N)-D(N) =:F(t,N), with N(0) = N m , 

where N, B,D € IR 7 . The numerical approximation of total number of particles in ith cell, 
Ni(t), is defined by iVj(i) which is the ith component of the vector N. The discretized birth 
term, Bf A , and death term, -Dp 4 , obtained from the cell average technique are defined below. 
These are the ith components of the vectors B and D respectively. Let us consider the Heaviside 
function 



( \ ifx>0, 



(2.13) H(x) :-- 

and 



if x = 0, 
if x < 0. 



(2.14) X: (x) - ' J '=' 



Xi Xi±i 

Then the birth and death terms are given as 

(2.15) Bf A : = B i ^ 1 Xr(v i _ 1 )H(v^ 1 -x l . 1 ) + B l X+(v t )H(v i -x l ) 

+BiX~{vi)H{xi - Vi) + B i+1 Xf (v i+ i)H(x i+ i -v i+ i), 
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and 

I 

(2.16) Df A := A = Ni(t)Y,P(?i,Xj)Nj(t)- 

i=i 

The first and the fourth terms on the right hand side of equation (|2.15p can be set to zero for 
i = 1 and i = I, respectively. The detailed formulation can be found in |12j . 
By using (|2,6p and (|2.7p the cell average technique (|2.1ip can be written as 

(2.17) = V^w-jw) 



< x j+ x k< x i-l/2 

+ [\f(v i )H(v i - Xl ) + \r(vi)H{xi - Vi)} 

X 

x i-l/2< x j+ x k< x i+l/2 

+\f(v i+1 )H(x i+ i -Vi+i) 



3>k , 1 v 

v> ■ _l_'y> . i-*" T . . _. ■ * * 



Xi+l/2 < x j + x k < x i+3/2 

-N^Y^xux^Njit). 
It should be pointed out here that in this work we consider the following discrete norm 

iiNii=xyi 



■I i 

8=1 



We consider C 2 ([o, b]) as a space of twice continuously differentiable functions on ]a, b[. Note 
that for the sake of simplicity in our analysis we assume that the coagulation kernel satisfies 

(2.18) e C 2 {]0,R]x]0,R}). 

Before moving to the main result, let us recall some definitions and a result taken from [1 lj . 

Definition 2.1. The local discretization error is defined by the residual left by substituting 
the exact solution N(i) into equation A2.12\) as 

(2.19) a{t) = ^ - (B (N(t)) - D (N(t))) . 
The scheme \2.12\) is called consistent of order p if, for Ax — > 0, 

\\a(t)\\ = 0(Ax p ), uniformly for all t,0<t< T. 

Definition 2.2. The global discretization error is defined by 

(2.20) e(t) = N(t)-N(t). 



6 



ANKIK KUMAR GIRI 



The scheme \2.12\) is called convergent of order p if, for Ax — > 0, 

||e(t)|| = 0(Ax p ), uniformly for allt,0<t< T. 

It is important that the solution obtained by CAT remains non-negative for all times. This 
can be easily shown by using the next well known theorem. In the following theorem we write 
M > for a vector M £ R ; if all of its components are non-negative. 

Theorem 2.3. Suppose that F(i,M) defined in \2.12\) is continuous and satisfies the Lipschitz 
condition as 

||F(t,P) - F(t,M)|j < L\\P - M|| for all P,MeR 7 . 

Then the solution of the semi-discrete system A2.12]) is non-negative if and only if for any vector 
M G W with M > ; the condition Mj = implies Fj(i,M) > for any i = 1, . . . ,1 and all 
t > 0. 

Proof. The proof can be found in [11, Chap. 1, Theorem 7.1]. □ 

Now we shall state the main result which helps us to show the convergence of CAT. Note that 
the proof of Theorem 12.41 is motivated by a convergence result in [T7] . 

Theorem 2.4. Let us assume that the Lipschitz conditions on B(N(t)) and D(N(t)) are satisfied 
for < t < T and for all N, N € M 7 where N and N are the projected exact and numerical 
solutions defined in \2. 5|) and \2.12\) respectively. Then a consistent discretization method is 
also convergent and the convergence is of the same order as the consistency. 

Proof. Using the equations (|2TT9D and (p^Ojl we have for e(t) = N(i) - N(t) 

j t e(t) = ait) + (B(N) - B(N)) - (D(N) - D(N)). 

We then take the norm on both sides to get 

j t Mt)\\ < \W(t)\\ + ||(B(N) - B(N))|| + ||(D(N) - D(N))||. 

Integrating with respect to t with e(0) = and using the Lipschitz conditions (|4.3p - (|4.4p we 
obtain the estimates 

< f || <t(t) II dr + 2L f \\e(r)\\dT. 
Jo Jo 



From this it follows by Gronwall's Lemma that 
(2.21) ||e(t)||<g[e X p(2Lt)-l], 
where 



eh = max \\o~ 

0<t<T 



If the scheme is consistent then lim/^o = 0. This completes the proof. □ 
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To fulfill the requirements of Theorem 12.41 f° r the convergence of CAT we need to prove 
that the scheme is consistent as well as the birth B(N(t)) and death D(N(i)) terms satisfy the 
Lipschitz conditions. 

3. Consistency 

In this section, the consistency of CAT for solving SCE (|1.2p is discussed. The discretization 
error is evaluated in the birth and death terms, respectively. Then the local discretization error 
is calculated by considering uniform and non-uniform smooth geometric grids. 

Let us begin with the integrated birth term of SCE (|1.2p over ith cell 

I rXi+1/2 r x 

Bi = - / (3(x -y,y)n(t,x -y)n(t,y)dydx. 

Z Jxi-l/2 JO 

By changing the order of integration we get 

1 ^ f x i + l/2 PXi+X/2 

Bi = tj2w / P(x-y,y)n(t,x-y)n(t,y)dxdy 



j = l"j-l/2 J x i- 1/2 
1 

2 



I r x i+i/2 r x i+i/2 
+- / / P(x -y,y)n(t,x - y)n(t,y)dxdy. 

«.-i/2 Jy 



1/2 

Now we apply the midpoint rule to the outer integrals in both terms on the right-hand side and 
use the relationship iVj = n(i,Xj)Axj + 0(Ax 3 ) for the midpoint rule to obtain 

1 

(3.1) B{ = — y jNjfy) / P(x — Xj,Xj)n(t,x — Xj)dx 

^ J Xi_i 10 



j = \ Jx i-l/2 
1 f x i+l/2 

+-Ni(t) / p(x-Xi,Xi)n(t,x-Xi)dx + 0(Ax 6 ), 

IXi 



2 



3\ 



=: Bi + 0(Ax 

Let us denote the integral terms in Bi by I\ and I2, respectively, and evaluate them separately. 

Integral term I\\ We consider the first integral term on the right-hand side in (|3.ip and use 
the substitution X X j — X to get 

(3.2) I 1 = -^Njit) / * 3 p(x',Xj)n(t,x')dx'. 

j = \ Jxi-l/2-Xj 

We now define lij and jij to be those indices such that the following hold 

(3.3) Si-1/2 - x j G Ai itj and 7^ := sgn[(xj_ 1/2 - Xj) - x h J, 
where 



sgn(x) := < 



1 


if 


X 


>0, 





if 


X 


= 0, 


-1 


if 


X 


< 0. 
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By the definition of the indices kj and jij in (|3.3p . the equation (|3.2p can be rewritten as 
(3.4) h = \^Nj{t) [ Xli < j+hi » ' /3(x',Xj)n(t,x')dx' 

i-1 k+lj + ^hi+l,}- 1 ) 



+\ E E T fc+1/2 ^ x >(*' ^o^' 

i=l fe=/ . . + i( Tii . + i) ■'"fc-i/a 



1 rXi+i/2-Xj 

1 ,-=1 . . 



J- 1 ! i+l,J + 3T*+l,J 

Let p be the total number of terms in the following sum 

ii+U'+KTi+lj- 1 ) rx k+1/2 

/ f3(x' , Xj)n(t, x')dx' . 

k=l i , j +Uj iJ +i) Xk ' 1/2 

In particular, let p := #{n : ijj + \(li,j + 1) < Ji < + oO^+lj ~~ -0} an d se ^ 

% := k,j + -^(lij + 1)j fc 2 := fci + 1, ■ ■ ■ , fcp-i := fci + (p-2). 

Next, we shall show that p is finite and can be estimated by a constant which is independent of 
the grid size. By using the definition of the indices lij and jij in (j3.3|) . we can estimate 

(p - 2)Ax min < Ax fc2 + Ax^ + . . . + Ax fcp _ 1 < ^(Axj + Ax i+1 ) < Ax 
which implies using the assumption of quasi uniformity (|2.ip that 



(p - 2) < Ax <i<:^p<i<r + 2. 

Ax m j n 

This means the above sum has uniformly bounded finite number of terms. So we can apply the 
midpoint rule to the integral in second term on the right hand side and use N^it) = n(t, x^)Axk + 
<D(Ax 3 ) to get 

"1 ^ — f X x 

(3.5) h = -E^'W / h ^' H " P(x',x 3 )n(t,x')dx' 

1 <_1 

+ 2^ N J^ E P(x k ,Xj)N k (t) 

J = 1 X i _i/ 2 <(x : j+X k )<X i+ i/2 

+-^N j {t) I 1 3 P(x',xj)n{t,x')dx' + <D(Ax 3 ). 

Integral term I2: Let us consider the second integral term in (|3.ip and use the substitution 
a; — Xi = x' to estimate 

I 2 = -Ni(t) P(x?,Xi)n(t,x?)dx. 
£ . n 
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Again by the definition of the indices kj and 7^ in (j3.3|) we split the above integral as 

^+1,1+5(72+1,1—1) 

k=l Jx k-l/2 

I'X 

o X.th I 



1 i i+1 ,i+2(7i+i,i-l) ~ Xk+1/2 

h = ^Ni(t) Yl / f3(x', Xi )n(t,x')dx' 

k=l Jx k-l/2 
I r x i+l/2~ x i 

+-Ni(t) / f3(x' ,Xi)n(t,x')dx' . 



S 'i+l,i+^7t+l, 



By applying the midpoint rule in the first term and using the definition of the indices Zjj and 

7ij, we get 

(3.6) I 2 = \Ni{t) P(x k ,Xi)N k (t) 

Xi+X k <X i+1 / 2 
1 r x i+l/2~ x i 

+ -Ni(t) / P(x',xi)n(t,x')dx' + <D(Ax 3 ). 



2 



By substituting (|3.5p . (|3.6[) into (|3.ip and using (|2.6p . we estimate 



Bi = 5 J] P(x k , Xj )N k (t) 

i =1 s*-i/2<( a, j'+ a, *)< a: *+i/2 

+9^*) S P(x k ,Xi)N k (t) 



2 

x i+ x k< x i+l/2 



p[x , Xj)n{t, x )ax 



-1/2~ X J 
fx;,-, /o-Xi 



1 1 x M 

+ 2 22 N i® I P(x', Xj )n{t,x')dx' 
+-A r 4 (t) / /3(x^Mi,£')^' + e>(Ax 3 ). 



1, . 1 



The first two terms on the right hand side can be combined as 

3>k 



( 1 \ 

Bi = Yl [l--5 j Ap{xkiX j )N 3 {t)N k {t) 

x i-l/2< x j+ x k< x i+l/2 

+-^JV 3 -(f) / P(x',Xj)n(t,x')dx' + 0(&x*). 



i— 1 ^+1,3+2^+1 
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Finally, by using (|2.6[) the above expression for Bi can be written as 



(3.7) Bi = Bi + ^Y^Njit) j h ' i+ ^ p( x ',Xj)n(t,x')dx' 

j = l x i-l/2~ x j 

+-ViVj(f) / P(x', Xj )n(t,x')dx' + 0(Ax 3 ). 

Let us denote the sum of the remaining two integrals on the right hand side in (|3.7p by the error 
Ei which will be discussed later. 

Now we concentrate to evaluate the integrated term Vi — XiBi by using ()2.3p and (|2.8|) as 
follows 

I rXi+i/2 r% 

Vi-XiBi = - / (x - Xi)/3(x - y,y)n(t,x - y)n(t,y)dydx. 

1 Jxi-i/2 Jo 

By changing the order of integration we get 

1 ^ f x J + l/2 f x i+l/2 

Vi-XiBi = / / (x - Xi )f3(x -y,y)n(t,x - y)n(t,y)dxdy 

I 

(x - Xi)f3(x - y, y)n(t, x - y)n(t, y)dxdy. 



j = l Jx j-l/2 Jx i-l/2 

\ r x i+i/2 r x i+i/2 



2 



x i-i/2 J y 



Now applying the midpoint rule to the outer integrals in both the terms on the right hand side 
and using the relationship JVj = n(t, Xi)Axi + 0(Ax 3 ) with (3(0, •) = 0, we obtain 

1 f^+l/2 

(3.8) Vi — XiBi = - Nj(t) / (x — Xi)j3(x — X j,Xj)n(t, x — Xj)dx 

2 j=l 

1 f x i+l/2 

+-Ni(t) / (x — Xi)f3(x — Xi,Xi)n(t,x — X{)dx + 0(Ax ), 



2 

=: Vi- Xi B t + 0{Ax b ). 

We denote the integral terms involving in Vi — XiBi by Pi and P2, respectively, and calculate 
them separately. 

Integral term Pi: Let us consider the first integral term in (|3.8p and insert to 
estimate 

(3.9) P 1 = ~Y J N j {t) / 1 ' {x' -Xi + Xj)P{x!,Xj)n(t,x')da!. 

j = l J^i-l/l-Xj 
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Py = J^Njit) / h - 3+r2ll ' ] (x' - Xi + x j )p{x',x j )n{t,x')dx' 

i-1 '»+i,j + |(7i+l,j-l) 

3=1 fc=/ 4 , Xfc " 1/2 

+->^A r j(t) / (V — Xj + Xj)j3(x , Xj)n(t, x')dx' . 



Zfc+l/2 

(ic' — Xj + Xj)f3(x' , Xj)n(t, x')dx' 



J— 1 ! i+lj'+TlTi+l,3 

Since the number of terms in the inner summation of second term on the right hand side is finite 
as before, therefore we can use the midpoint rule to the integral in second term on the right 
hand side and use iVfc(t) = n(t, x^Ax^ + 0(Ax 3 ) to obtain 

-j 1 1 /*X ]_ 

(3.10) Pi = ^^2Nj(t) / h ' 1+ ^ 1 ' 3 {x' - Xi + x^^x' , Xj )n{t,x')dx' 

j=l Jxi-l/2-%j 



J 

i-1 



1 

3=1 a:i-i/2<fe+a;fc)<a; l +i/2 



i-1 



1 fXi+l/2-Xj 

+-2^iVj(t) / (rf -Xi + Xj)P{a!,Xj)n(t,rf)dx' 

3=1 X Ij+ij+^Ti+l j 

^ i-1 ^+l,j + |(7i+l,j-l) _ 

+ 24 Z)^) £ Ax fe 3 — {/?(x fc ,x>(t,x fe )} 

i =1 fc=^,j+^(7*,j+i) 
+C(Ax 4 ). 

Integral term i-*2 : Let us consider the second integral term in (|3.8p and use the substitution 
x — Xi = x' to estimate 



2 2 



I f x i+l/2~ x i 

-Ni(t) / xfp(a/ f Xi)n(t,a/)dx. 
1 Jo 



By the definition of the indices and 7ij in ()3.3p we split the above integral as 

Ji+l.i+KTi+l.t-l) p fc+1/2 
P 2 = rJV<(t) £ / x'p(x',Xi)n(t,x')dx' 

k=l Jx h-l/2 
]_ r x i+l/2~ x i 

+-Ni{t) / x'f3(x',Xi)n(t,x')dx'. 

'i + l,i + 7">i + l,i 
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We apply the midpoint rule in the first term on the right hand side to obtain 



(3.11) P 2 = -Ni(t) (xi + x k -Xi)P(x k ,Xi)N k (t) 

Xi+X k <X i+l / 2 



1 f x i+l/2~ x i 

--Ni{t) I x'P(x',Xi)n(t,x')dx' + 0(Ax 4 } 



! <+l,<+2Ti+l,j 

By substituting (|3.10p . and (|3.1ip into (|3.8p . we estimate 



Vi-XiBi = -^2Nj(t) ^2 (xj + x k - Xi)P(x k ,Xj)N k (t) 

J = 1 X i _ 1/2 <(Xj+X k )<X i + 1 /2 

+ ~ N i(t) XI (Xi + X k ~ Xi)l3(x k ,Xi)N k (t) 



2 

Xi+X k <X i + 1/2 

I x, . . , 1 . 



j—l J x i — ±/2~ x j 

+^2^Nj{t) / {rf -Xi + Xj)p{rf,Xj)n(t>a!)da! 



J l i+l,j+5T»+l,i 

i-1 'i+i,i+a(Ti+M _1 ) 



+ ^X Ar i( t ) E Ax fc 3 — {/?(x fc ,x>(t,x fc )} + 0(Ax 4 ) 

J =1 fc=ii,i+|(7»,j+l) 



The first two terms on the right-hand side can be combined and the above equation can be 
written as 



Vi-XiBi = ^ (l--5 jik j(xj+x k -x i )(3(x k ,x j )N j (t)N k (t) 

Xi-l/2<Xj+X k <X i+1/2 

j = l J x i—l/2~ x j 

+-2_^Nj(t) / (x' - x i + Xj)/3(x , ,Xj)n(t,x')dx' 

2 A, , i 



J l i+l,j+!fT»+l,i 
i-1 'i+l,j + 5(7i+i,j-l) 



+ 24 Z)^*) X Ax fc 3 — {/3(x fc ,x>(t,x fe )} + 0(Ax 4 ). 

•? =1 fc=^,j + |(7i,j+l) 
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Finally, by using (|2.9[) and (|2.6p . the above expression can be rewritten as 
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(3.12) Vi-XiBi = Vi-XiB, 

i-l 



1 f 

.7=1 JX 



X l- ■+ 1 7- 

l ' J l,J (x — Xi + Xj)f3(x', Xj)n(t, x')dx' 

i-\/2~ x j 



1 * fXi+l/2-Xj 

-^iVj(t) / (x' - Xi + x j )(3(x',x j )n(t,x')dx' 



i=i 



a;, , l 



^ i-l ^+i,i + |(7i+i,j- 1 ) q 

+ 24 Z^*) S Ax fc 3 — {/3(x fc ,x>(i,x fc )} 

+C(Ax 4 ). 



Now we evaluate the each term in (|2.15p separately. We begin with the first term without 
Heaviside function H{x) and insert the value of A from (|2.14p to get 



\~(v- {)B- i = Vi ~ 1 ~ Xi ~ l B i = 

Xi- Xi-i 1 Axi + Axj_i 



Using the equation (|3.8|) and (|3.12p . we obtain 



(3.13) Xrfa-JBi-t 



Axi + Axj_i 

i-2 



Vi-i - Xi-xB, 



i-l 



-y^iVj(t) / 1 1 ' J+ ^"' 1 1,3 (x' -Xj_i + Xj)(3(x',Xj)n(t,x')dx' 



i-l 



2_^Nj(t) I (x - Xi-x + Xj)P(x ,Xj)n(t,x )dx 



3=1 
i-2 



■sE^w E 

j=i fc=i i _i ij +|( 7i _i J +i) 



+C(Ax 4 



Ax fc 3 — {/3(x fc , Xj)n(t, x fc )} 



To solve equation (|3.13p . we estimate V^_i — Xj_i-Bj_i by using (|3.8p as follows 



Vi-i - Xi-xB. 



i-l 



i-2 



i=i 



-1/2 



(x — Xj_i)/3(x — Xj,Xj)n(t,x — Xj)dx 



-1/2 



+ I (x - Xi_i)/3(x - Xj_i,Xj_i)n(t,x - Xj„i)n(t,Xj„i)Axj_i(ix]. 
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By setting f(-,y) := f3(-,y)n(t, •), the above equation becomes 

i-2 



(3.14) Vl-i-Xi-iBi- 



-1/2 



(x — Xi-x)f(x — Xj,Xj)dx 



-3/2 



(x - Xi-x)f(x - Xi-i,Xi-i)n(t,Xi-i)Axi-idx\. 



We use Taylor series expansions of each integrand about Xj_i in equation (|3.14p as 

(x - Xi-l)f(x - Xj,Xj) = + /(xj_l - Xj,Xj)(x-Xi-l) 

+f x (x i - 1 - Xj,Xj)(x - Xj_i) 2 + 0(Ax 3 ), 
(x - x;_i)/(x - Xi-i,Xi-x) = + /(xj-i - x i -i,x i ^ 1 )(x - Xj_i) + 0(Ax 2 ). 

The substitution of the above Taylor series expansion in equation (|3.14p gives 

1 T 1 ^ 2 

~ T^^Z N i fx ( x i-i~ x 3 > x i ) Aa; i-i 



Vi-i - Xj_ii?i_i 



1 



+ o/( x '*-i - Xi_i)n(t, Xi_i)Ax i _ 1 + 0(Ax ) 

o 

Since /3(xj_i — Xj_i,Xj_i) = /3(0, Xj_i) = 0, therefore we have /(xj_i — Xj_i,Xj_i) 
implies that 

i-2 



0. This 



Vi-i - x^ 



1 

= 24 E^^- 1 " x i> x j) Ax ti + G(Ax 4 ). 



Again the application of Taylor series expansion gives us 

1 i_2 

(3.15) V5_i - Xi-ifli-i = -^JVj/^ij - Xj-.x^Ax?.! + 0(Ax 4 



Finally, substituting (|3.15|) into (|3.13p . we obtain 



1 i_2 

(3.16) A- (vi-i)Bi-i = — Njf x (xj - Xj,Xj) 

i-2 



Axf_ x 



Axi + Axj_i 

' ( i— l,j + 2Ti— i,j (x' — Xj_x + Xj) 



i=i 

i-2 ^,j + |(7i,j-l) 

^5>w E 



f(x',Xj)dx' 



x ! Axi + Axj_i 



Axj + Axj_i 

f(x',Xj)dx' 



Ax fe 3 



12^" JV " / ^ Axj + Axi_i 

J=l fc=2 i _i ij +§(7 i _i i j+l) 

+0(Ax 3 ). 



/a;' > -Ej ) 
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Next, the second term in (|2.15p is evaluated as 
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Ax i+ i + Axi 



Calculating as before, we estimate the above expression in the following form 



(3.17) 



1 <_1 

Xf(vi)Bi = Bi- — ^ Njf x (xi - Xj,Xj) 



Ax? 



Axi + Ax i+ i 



i-i 



3=1 



Axi + Ax i+ i 



» r / s f Xi+1 / 2 ~ Xj ix' -Xi+Xi) r , 

+ £#i( t ) I — , A „ 31 f(x', Xj )dx' 



3=1 



Axi + Ax i+1 



1 i-1 ii+XJ+|(Ti+Xj-l) ^ 3 

J=l fc=«ij4-|(7i,j+l) 

+0(Ax 3 ). 



Similar to the second term we obtain 



(3.18) 



1 l ~ l 

\ {vi)Bi = Bi + — ^ Njf x (xi - Xj,Xj) 



Ax ? ; 



3=1 



Axi + Axi-i 



i-l 



3=1 



^itj^^^iij (^ X% X 



1 ' 3 1 f ( I 



Axi + Axi-i 



f(x', Xj)dx' 



ST- v f x i+l/2~ x 3 (x 1 - Xi + Xj) , , 



3=1 



x , i i 

'i+X,3 + 27i+l, 3 ' 



Axj 4- Axj_i 



i-l ii+xj+|(Ti+xj-l) a 3 

-±y N f t ) y — — 

+C(Ax 3 ). 
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Finally, similar to the first term we can easily estimate 



1 - Ax 3 
(3.19) \f(v i+1 )B i+1 = — N ifx{xi - Xj,Xj) Ax _ - l ~£ x . +i 



j= l J Xi+l/2-Xj 

» T / n r^ 2 ~ X] (X' ~ X l+1 + Xj) r . , . 1 . 



3— 1 X 'i+2,j + h~<i+2 



3 

h+2,j + \{li+2,j-l) 



+^X>0 E Ax^A^/*'^'^ 

3 k=k+l,j+2 (7i+l,j+l) 

+0(Ax 3 ). 



By substituting (fXTUj) . ([3TTT1) . (pTTgj) and (IXlH into (l2TToD and using ([3771), the local discretiza- 
tion error can be evaluated as follows 

Case I: Vi-i > xi—i, Vi > Xi and Vi+i > Xj+i 



(3.20) 



+ 



1 



5 



CA 



5, 



j = l J^z-l/2-^j 



Ax 3 _ x 



• r*i+i/2-*j 
Ax? 



/(x', Xj)dx' 



12 \ Ax,, + Axj_i Axi + Axj+i 
i-i 



3=1 Jx ' 



E^ 



x . Axi + Ax i+ i 



f(x', Xj)dx' 



V- « r , v / ,<B *+ 1 /»-*i (x' - ^ + x 7 -) , . , 



3=1 



i-2 



Axi + Axj+i 



E, r /,\ f Xl i-l,3 + hli-l,3 i x ' ~ x i-l+ x j) t . I , i I 
Njit) I —. r —j(x,Xj)dx 

Jx ,,n-x- Axi + Axi_i v v 

j = l Jx i-3,/2 x ] ' 



t-1 



E„r / N Z^ 1 " 1 / 2 ^ (X' - Xi_l + Xi) „ , , , 

N jit) I a_ , IT 31 f{x', Xj )dx' 



3=1 



x i- 



Axj + Axi_i 



So 
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i-l 



4e« e 

i =1 fe=ii,j + 5(7i,j+l) 

sS« E 



Axj + Ax m 
Ax fc 3 



fc=/f_lj + i(7 i _l ] j+l) 



+0(Ax 3 ). 



fx'(%ki Xj) 
fx' { x ki Xj , 



> =: E 4 



Case II: Vi-i < Xj_i, Uj < Xj and < Xj+i 
Similar to the previous case, we have 



(3.21) 



D 



CA 



i-l 



3=1 



x i+l/2- x j 



+ 



'2^2^ J x f(x', Xj )dx' 

3 = 1 ij + lj + 27i+lj 

Axf +1 \ 
12 VA^ + Ax^i Axi + Axi+J^ 



Ax? 



Njf x '{xi — Xj, Xj) > — : E' 2 



+1 ' j+hi+1 » {x '- Xi+1+Xj) f(x>,x 3 )dx' 



+3/2-XJ ( x ' _ x + x ) 

j(x , Xjjax 



i+2j + 2Ti+2J 



Axj + Axj + i 



+E^w / 

i— 1 „ 

-E^-w / 
-E>*(*) 



l i,j + \~<i,j (X Xj + Xo) . / \j / 

/(x , Xjjax 



_ x Axi + Axj_i 



'i+l/2-Xj (x'-X i +X j ) , , 

j(x,Xj)dx 



i+2,j + §(7i+2,j-l) 



4J>'"> E 

i=l fe=Z i+ i, :7 - + i(7 i+ i J -+l) 

i-l '»+i,j+5(7»+i,j-1) 

E 



Axj + Axj_i 
Ax fc 3 



Axj + Axj + i 
Ax fc 3 



fx 1 {Xki Xj) 



3=1 



k=l i,j + h{li,j+l) 



Axi + Axj_i 
+C(Ax 3 ). 



fx' (Xk i Xj 
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Case III: v j_i < Xi-x, V{ = Xi and Vi + \ > Xi + \ 
(3.22) 



B? A = * 



i-l 



/(a/, Xj)dx' 



a*- 1/2 —ay 



=: Si 



"aE^iW _/ f(x', Xj )dx' 

i = l 'i+lj+^Ti+l.i 

+C(Ax 3 ). 

Next, the discretization error for death term is calculated in the ith cell. From equation ([2? 
the integrated death term can be written as follows 

Di= } j l K(x,y)n(t,y)n(t,x)dydx. 

^ x i-l/2 j = l ^ x j-l/2 

The application of the midpoint rule to the outer and inner integrals gives us 

I 

(3.23) A = Ni(t) ^ K(xi, Xj)Nj(t) + 0(Ax 3 ) = A + 0{Ax 3 ). 

i=i 

From the equations (|3.20H3.22j) and (|3.23p . we can estimate the local discretization error &i{t) 
(Bi - A) - {B? A - Bf A ) as 

'E 1 + E 2 + E 3 + E 4 + 0(Ax 3 ) ifi€Oi, 
E 1 + E' 2 + E' 3 + E' 4 + 0(Ax 3 ) ifi€a 2 , 
A + C(Ax 3 ) if i G 2t 3 . 



(3.24) 
where 



2d ={iG N|«i_i > Xi-i,Vi > Xi,v i+ i > Xi+i}, 
2l 2 ={iG N|«i_i < Xj_i,Wj < Xi,Vi+i < x i+ i}, 
% = {i € N | < Xi-i,Vi = Xi,v i+ i > x i+ i}. 

Here we consider three different cases to find the order of consistency. Then, the order of 
consistency is given by 

(3.25) ||a(t)|| = + E + E 

ie% x iesu ie2t 3 

The following types of grids will be considered to find the order of consistency of CAT. 

3.1. Uniform grids. Let us begin with the case of uniform grids i.e. Axj = Ax and Xi = 
(i — 1/2) Ax for any i = 1, ... , I. Here, Ei and E' 2 defined in (|3.20p and (|3.2ip . respectively, 
are obviously zero. In case of such uniform grids, we have 

x i~l/2 ~ x j = x i—ji x i+l/2 ~ x j = x i—j+li an d x i-3/2 ~ x j = x i— j— 1- 

By using the definition of indices Zjj and 7^ in (|3.3p . we calculate 



c i— 1/2 
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which gives 

x i—l/2 ~ x j = x i—j = x h,j' 

Similarly, we obtain 

x i+l/2 ~ x j = x i—j+l = x l i+ iji 

and 

x i-3/2 ~ x j = x i-j-l = x k-i,j- 

This shows that 7i-ij = Ji t j = 7i+ij = 0. Therefore, in (|3.20p - (|3.2ip . the error terms E±, 
E% and E' 3 become zero. It can also be easily realized that Xj_ 3 / 2 — Xj, Xj_iy 2 — Xj, and 
x i+i/2 ~ x j are the pivot points of the adjacent cells, i.e. k-ij = (i — j — l)th, kj = (i — j)th and 
k+ij = (i — j + l)th cells, respectively. Thus, by substituting the values of all these indices in 
E4 and £4 defined in (|3.2U|) and (|3.2ip . respectively, and applying the Taylor series expansion, 
we obtain E 4 = 0(Ax 3 ) and E' 4 = 0(Ax 3 ). Then, from IpTM)) . we obtain 

<7i(t) = 0(Ax 3 ) if i G 2t x ,2l 2 ,2l3. 

By using ()3.25j) , the order of consistency is thus given by 

\\a(t)\\ = 0(Ax 2 ). 

Therefore, the cell average technique is second order consistent on uniform grids. 

3.2. Non-uniform smooth grids. Non-uniform smooth grids can be obtained by applying 
some smooth transformation to uniform grids. Assume a variable £ with uniform grids and a 
smooth transformation x = g(£) such that x i± i = for any i = 1,...,J to get non- 

uniform smooth grids. In this case, we show that the scheme is first order consistent. Let 
h be the uniform mesh width in the variable £. For such type of smooth grids, Taylor series 
expansions in smooth transformations give 

Axi = x i+ x - x t _i = g{ii + ^) - g{ii - ^) = hg'(&) + 0(h 3 ). 
Hence, by calculating Axj_i and Ax,+i similar to Axj, we obtain 

Axi - Axi-i = 0(h 2 ), 

A Xl + Ax i+1 = h[g'(Ci) + </(&+l)] + 0{h 3 ) = 2hg'(Ci) + 0(h 2 ), 
and similarly, we have 

Ax, + Ax^ = h[g'(Zi) + ff'fe-i)] + 0(h 3 ) = 2hg'(Zi) + 0(h 2 ). 

In particular, we deal with a special type of non-uniform smooth grids which is known as 
geometric grids. Such type of grids can be defined as x i+ i = rx^i, r > 1, i = 1, . . . ,1. An 
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exponential function can be applied on uniform grids as a smooth transformation to construct 
such type of geometric grids. Mathematically, we write 

2*4-1/2 = exp(£ m/2 ) = exp(/i + Ci-1/2) 

= exp(/i)exp(^_ 1/2 ) 

= exp(/i)x i _ 1/2 =: rXj_ 1/2 , r > 1. 

To solve the error terms appearing in (|3.24|> . let us further assume that £n, £i 2 , £21, £22, £31 and 
^3 2 are corresponding points on uniform mesh for x t . +i j . +1 + i 7 . +1 2^+1/2 — x j+i, x^.^i^.., 
x i-i/2 ~ x ji x i i _ 1 3 _ 1 +i 7l _ 1 j _ 1 an d a^i-3/2 — x j-ii respectively. Due to an application of expo- 
nential smooth transformation, these points can be defined as 

in = ln(^. +li . +1+ i 7i+1 . +1 ), . . . ,£32 = rn(^_ 3/2 - Xj - X ). 

By the definition of the indices in (j3.3[l . we know 

Xi+1/2 ~ Xj+l e A Zi+lj+1 , ^i-l/2 - »j G Aj^ and Xi _ 3/2 - G A^^. 

For geometric grids, we have 

Xi+l/2 ~ Xj+l = r{xi_ 1/2 - xj) = r 2 (xi_ 3/2 - Xj-i). 

Therefore, we have 

^i+lj'+l = h,j + 1 = h-i,j-i + 2. 
Further, in case of geometric grids, we have 

7i+l,i+i = = Tt-iJ-l- 

Let us consider 

hi = £11 - £12 = ln(x z . +i . +i+ i 7 . +ii . +i ) - ln(s i+1/2 - x j+ i) 



2 

x i+l/2 — / V^i— 1/2 

^i-3/2 - SCj-l / 



Similarly, we estimate 



and 



ei2-62 = ln( Xm/2 " Xj+1 )=ln(r) = /., 

x i-l/2 ~ x j 



62 - 62 = ln( Xi 1/2 Xi ) = ln(r) = h. 

\ x i-3/2 — x j-l J 

Again, by application of smooth transformation, we can easily obtain 

(3.26) x h+h , +1+ i 7i+lj , +1 - (x i+1/2 - x j+1 ) = 5(61) - 5(62) = M'(6 2 ) + 0(h 2 ), 

(3.27) . " fo-i/2 " Xj) = gfai) ~ 5(62) = + ©(/j 2 ), 
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and 

(3.28) ^_ i ._ i+ i 7i _ ij ._ i - - xj-i) = 5(61) - 0(62) = W(&2) + 0(/i 2 )- 

All these identities will play an important role to solve the error terms involved in (|3.24p . which 
helps us to calculate the order of local discretization error o-j. 
We first evaluate E\ as follows 

Applying the left and right rectangle rules in the integrals involved in the first and second terms, 
respectively, on the right-hand side, we estimate 



1 

i=i 
1 - 

+^^'(^^1/2 ~ ^'>^)(^+i/2 - ay - ^ i+lj . + i 7i+lj .) + 0(Ax 2 ). 

Then an application of Taylor's series expansion about a^-i/a = £1+1/2 in f{x i _i/ 2 ~ x j) gives 

1 £i 

i=i 
1 - 

-jE^'O/^ ~ x 3> x 3)( x i i+1J +i<y i+1J ~ Xi +V2 + x i) + 0(Ax 2 ). 
i=i 

We replace j by j + 1 in the second term on the right-hand side and use the relationship 
Nj(t) = n(t, Xj)Axj + C(Ax 3 ) for the midpoint rule. Also, we drop the term which is of second 
order, and obtain 

1 i_1 

Ex = - X j) Ax jf(%i+l/2 - Xj, Xj)^.^!^, - Xi _ 1/2 + Xj) 

3=1 



1 

- ^ n(t, x j+1 )Ax j+1 f(x i+1/2 - x j+1 ,x j+1 )(x k+i . +i+ i 7i+i . +i - x i+1/2 + x j+l ) 
3=1 
+0(Ax 2 ) 



Approximating the function x h-> n(t, x)f(xi±i/ 2 — x, x) at Xj by n(t, x)f(xi±i/ 2 — x, x) evaluated 
at x = Xj + \ in the first term, we evaluate 



1 

(3.29)Ei = - Y. \ A - r J { - r > ~ X *-V2 + x j) ~ Ax 3+l( x k +hj+1 +± Ji+ld+1 ~ x i+l/2 + Xj+l)} 

3=1 

xn(t,Xj+i)f(x i+ i/ 2 - x j+1 ,x j+ i) + 0{Ax 2 ). 
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By using the identities in the beginning of this section, we calculate 

(3.30) Ax, (xi itj+ l*Yij ~ X i-V2 + x i) ~ Ax J+i( x k +1 , J+1 +^ +1 , J+1 ~ x i+i/2 + x j+ i) 
(Axj - Ai j+ i)(x Iij+ i 7ij - Xi _ l/2 + Xj ) 

-Ax J + 1 [(x h+i j + i + hi+i j + i - X t+1/2 + X j+1 ) - {X kj+ l 7ij ~ Xi-l/2 + Xj) 

£7- 



0(h 2 ){h ig '^ 22 ) + 0(h 2 )} - {hg'{i j+l ) + Oih^jilnWitu) - </(£ 22 )} + 0{h 2 )] 



0(h 3 ). 

Therefore, substituting ()3.30p in (|3.29p . we obtain 
(3.31) E x = 0(Ax 2 ). 

Next, let us calculate E% defined in (|3.20p as follows 

-3 = tmQ^^^,^ 



+ yNAt) / \ ~r — OJ -f[x,Xj)dx 

i=! Jx, ^ A Xi + Ax l+1 Jy 3 > 

YNJt) / h ~^ + ^ {X - Xi -1 + X i> f( x , xj )dx 



j = l ->^i-S/2-Xj 

- > iV,-(t) / — —f(x,Xj)dx. 

t< 3 Jx, ^ Axi + Axi-r J > 

Applying the left rectangle rule to the integrals appearing in first and third terms, and the right 
rectangle rule to the integrals in second and fourth terms, we estimate 

1 il Ax- 

E 3 = ~2 Y, N dt)f( X i-l/2 ~ Xj, Xj) ^ _ fa. . + hi> . - Xi _ 1/2 + Xj 

i=i 



-\^2 N j( t )f(x l+ l/2 -Xj,Xj)^-^—(x h+ij + l^ +h3 -X i+1/2 + X j) 

i=i ' 1 

1 Ax- i 

+ 2 E ^^)/(^-3/2 ~ gj, Sj) A + t A ^_ i (^ i _ 1 ,,+| 7i _ 1 , j ~ gi-3/2 + Xj) 
j=l * 1 

1 ' 1 Ax- 
+ 2E JV iW/( 3; w/2 - x i' x i) Ax . i^ t^+hj ~ x *-i/2 + x i) + 0{Ax 2 ). 
j=i 



Let us approximate / at {x^^i 2 — Xj, Xj) by / expanded around (#j-i/2 — Xj,Xj) in the third term 
and / at {xi-\j 2 — Xj,Xj) by / expanded around (x i+1 / 2 — x j,Xj) in the fourth term. Further, 
we replace j by j + 1 and j — 1 respectively in second and third terms. Also, the relationship 
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Nj(t) = n(t,Xj)Axj + 0(Ax 3 ) is used to get 

1 %1 Ax- 
E 3 = x 3 )Ax 3 f{x^ 1/2 - Xj ,Xj) A ^ - ^— O^.+f 7iJ ~ Xi-i/2 + gj) 



X- 



1 

--J^n(t,a;j + i)Axj + i/(x i+1/ 2 - x j+ i,x j+ i) 
i=o 

AXj . . 

Ax i + Ax m ° E ^+ 1+ ^+^+ 1 " Xm/2 + Xj+l) 
1 i_1 

+ 2 X] ^-1)^-1/(^-1/2 - Zj-i, ^i-i) 

J'=2 

Axi_i 

X-: : X, ,1 , — Xjv/o + Xj-i ) 

Axi + Axi-i 'i-l,j-l + 27i-l,i-l « j 17 

1 *~ 2 Axj ] 

+2 5^«(t.^) A «j7(»i+i/2 - x j' x j) ax~+^~7 ^+2^ ~ Xi ' 1 ' 2 + Xj ^ 

j=l 1 1 

+0(Ax 2 ). 

Without loss of generality, we can drop the terms which are second order accurate. Moreover, 
we approximate the functions x i-> n{t,x)f{x iT i/ 2 — x,x) at point Xj by n(t, x)/(xj Tl/ / 2 — x,x) 
evaluated at points x = Xj T \ of the first and fourth terms, respectively, to obtain 

j=2 1 1+1 



-^2n(t,x j+1 )f(x i+1/2 - x j+1 ,x j+ i) 

3=1 

AxjAxj+i 
Axj + Axj+i 

1 

3=2 



x • +A K + i, 3+ i+h <+ i, J+ i " X *+V2 + 



Axj_iAxj_i 

^ (X; ,1„ — X,_q/o + Xj-l) 

Axi + Ax,_i v ii-i,i-i+27i-i,i-i * •v 2 ^ •? ^ 



X- 

1 ** A> 

3=1 

„2^ 



+0(Ax z ). 

Let us denote each summation with the factor \ on the right-hand side by En, . . . , Eu respec- 
tively. Therefore, we can write 

(3.32) E 3 = (£33 - #3i) + (£34 - £32) + 0(Ax 2 ). 
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To simplify (|3.32p , we first calculate E 33 — £31 as follows 



1 

(3.33) £33-^31 = -^2n{t,x j -i)f(x i _ 1/ 2 -Xj-uxj-.!) 



i=2 

X 



Axi-iAx,--! 

-; 7—* — (x, , i„, — Xavi? + Xj-i) 

Axi + Axi-i h-ij-i+iii-ij-i * 1; 

Ax{Axj 

< x k. i +hi, i - x i-vz + x ^ 



Axi + Axj+i V '*.J + 3Ti,i 

Again by using the identities mentioned in the beginning of this section, we need to estimate 
the following term for solving f|3.33j) . 

Axi-xAxj-iiAxi + Aa; i+ i)(x l( _ liJ _ 1+ i 71 _ lii _ 1 - Xi_ 3/2 + a^-i) 
-AxiAxj(Axi + Aa; i _i)(x IiiJ+ i 7ljj - x;_ 1/2 + Xj) 
= Axj_i (Axj_i - Axj) (Axj + Ai i+1 )(s,._ i) ,_ i+ i 7 ,_ i) ._ i - + xj-i) 

=o(h 2 ) 

+Ax j [Ax i - 1 (Ax i + Ax,+i)(x /i _ i j _ i + i 7i _ ij _ i - Xl _ 3/2 + Xj_i) 
-Axi(Aa?i + Ax i _i)(x,._ i+ i 7 ._ . - x;_ 1/2 + 

= o(/i 5 ) + (/ i5 '(e J ) + o(/ i 3 ))[( %'te-i) +o(/ l 3 ))(2/ l5 'te) + o(/ l 2 ))(M'fe 2 ) + o(/ i 2 )) 

=V(Ci)+c(/i 2 ) 

-(%'(&) + 0(h 3 ))(2hg'(^) + 0(h 2 ))(h ig '(t 22 ) + 0{h 2 ))} 
= OQC)[2h 2 {c/{^) 2 h x y{^) - </(£ 22 )}] + 0(h 5 ) 
= 0{h) ■ 2h 3 h 1 (g'(^)) 2 g'^ 32 ) + 0(h 5 ) 
= 0(h 5 ). 

Inserting this estimate in (|3.33p . we obtain 

(3.34) E 33 -E 31 = 0(h 2 ). 
Analogous to (|3.34p . we can easily show that 

(3.35) E 3A -E 32 = 0(h 2 ). 
Finally, substituting (|3.34p and ()3.35p into (|3.32p . we have 

(3.36) E 3 = 0{Ax 2 ). 
In a similar way, we can prove that 

(3.37) E' 3 = 0{Ax 2 ). 

Next, it can be easily observed from (|3.20p and (|3,2ip that the error terms E 2 , E' 2 , E A and E' 4 
are second order accurate independent of meshes. Therefore, by substituting (|3.3ip . (|3.36p and 
(l3T3Tjl into (pM) . we have 

Gi(t) = 0{Ax 2 ) if % G %i,&. z ,%. 
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Thus, using (|3.25p . we obtain 

\\*(t)\\=0(Ax). 

This shows that the cell average technique is first order consistent on such type of non-uniform 
smooth grids. 

Remark 3.1. It should be pointed out that, due to the cancellation of second order terms, the 
error terms E2, E' 2 , E3, E' 3 , E4 and E^ can be shown third order accurate on geometric grids. 
However, since this will not improve the order of consistency (because E\ is only second order 
accurate for such grids), we do not include further calculations. 

4. LlPSCHITZ CONDITIONS ON B(N(i)) AND D(N(t)) 
Let us consider the birth term for < t < T and for all N, N € R J . We get from (|2,15p 



|B(N)-B(N)|| < J2K^i-i) H ^i-i- x i-i)\Bi-l(^)-Bi-x(N)\ 

i=l 

I 



i=i 

+ ^A+(TJi + i)F(^ +1 -u i+1 )|S i+1 (N)-B m (N)|. 
i=i 

The definitions of \ (x) and H(x) in (|2.14p and (|2.13p . respectively, guarantee that < 
\f{x)H(x) < 1. Thus, by using this upper bound, the above inequality becomes 

(4.1) ||B(N)-B(N)|| < ^2 l-^i-l(N) — J3i_i(N)| + ^ l-Bj(N) — -Bj(N)| 

i=\ i=l 
I 

+ £|J3 i+ i(N)-fl i+ i(N)|. 

i=i 

By (|2.18p . there exists a constant C > such that /3(x,y) < C for all x,y g]0, R]. Then, 
substituting (|2.6p into (j4. ljl . we have 

I i-l 

||B(N)-B(N)|| < ^EE E \NjNk-NjNkl 

i=l j=l X i -3/2<Xj+X k <X i _ 1 / 2 
1 1 * 

i=l j=l X i ^ 1/2 <X J +X k <X l+1 /2 

1 I i+l 

+2 C EE E i^-^Ai 

8=1 j=l x i+1/2 <a: j +x fc <a: i+ 3/2 

3 ' 7 

^ ^EEi^-^i- 

i=i fe=i 
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Now we apply the following useful equality NjN k - NjN k = |[(iVj + Nj)(N k - N k ) + (Nj - 
Nj)(N k + N k )] to get 

3 1 1 r 1 
(4.2) ||B(N)-B(N)||<-CX;E 1(^ + ^)11(^-^)1 + 1(^-^)11(^ + ^)1 • 

j=l k=l L J 
It can be easily shown that the total number of particles decreases in a coagulation process, i.e. 
I 

Nj < iVj. := Total number of particles which are taken initially. 

i=i 

The equation (|4.2p can be rewritten as 

■ I I 

Y,m-N k )\ + Y,\(Nj-N 3 

k=l j=l 

< 3iV£(7||N-N||. 



(4.3) ||B(N) -B(N)|| < -N T C 



Similarly as before we can easily show the Lipschitz condition for death term as 
(4.4) ||D(N) -D(N)|| < 3i\#C||N-N||. 

Thus, Theorem 12.41 implies the convergence of the cell average technique and the convergence is 
of the same order as the consistency. 

5. Conclusions 

We have presented a detailed convergence analysis of the cell average technique for nonlinear 
continuous Smoluchowski coagulation equation. It is proved that the cell average technique is 
second order convergent on uniform grids. However, it gives only a first order convergence on 
non-uniform smooth geometric grids. To obtain a second order convergence, either one needs 
to adapt a different approach than the one presented here, or modify the error term E\ which 
may lead to some improvements in CAT. It is also interesting to analyze CAT for nonlinear 
continuous SCE on more general grids, which we intend to study in future. 
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